Polarized Cosmological Gravitational Waves from Primordial Helical Turbulence 
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We show that helical turbulence produced during a first-order phase transition generates circularly 
polarized cosmological gravitational waves (GWs). The characteristic frequency of these GWs for 
an extreme case of the phase transition model is around 10 -3 — 10 -2 Hz with an energy density 
parameter as high as 10 -12 — 10 -11 . The possibility of detection is briefly discussed. 
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Since cosmological GWs propagate without significant 
interaction after they are produced, once detected they 
should provide a powerful tool for studying the early Uni- 
verse at the time of GW generation [lj. Various mech- 
anisms for cosmological GW generation have been stud- 
ied, including: quantum fluctuations during inflation Q; 
bubble wall motion and collisions during phase transi- 
tions 0; cosmological magnetic fields 0,13; and plasma 
turbulence 0,0,0- 

In this letter we focus on polarization of cosmological 
GWs generated by helical stochastic turbulent motions 
[E IS 03 • We find that helical turbulence generates cir- 
cularly polarized stochastic GWs and we compute the 
polarization degree. The formalism we use is general and 
can be applied to study the generation of stochastic GWs 
by any helical vector field (e.g., helical magnetic fields 
PJJ,[ll|]). Primordial polarized GWs might be generated 
from quantum fluctuations accounting for the gravita- 
tional Chern-Simons term |l2j . 

GWs are sourced by the transverse and traceless part 
of the stress-energy tensor T^ v In the case at hand 
T^ v describes a turbulent cosmological fluid after a phase 
transition 0. 0- For spatial indices i ^ j, I\y(x) = 
(p+p)iti(x)uj-(x), where p and p are the fluid pressure and 
energy density and u(x) is the fluid velocity. The fluid 
enthalpy density p + p is taken to be constant throughout 
space. The transverse and traceless part of in Fourier 
space is Iltf(k) = [Pu(k)P jm (k) 



where P:,(k) = Sn—kikn with ki 



iPy(k)P iro (k)]T ;ro (k), 

-ki/k. Consistent with 



observations we have assumed flat space sections. 

To model the turbulence we assume that in the early 
Universe at time t- m (at a phase transition) liberated vac- 
uum energy p vac is converted into (turbulent) kinetic en- 
ergy of the cosmological plasma with an efficiency k over 
a time scale T s ti r on a characteristic source length scale 
Ls Q. After generation, the turbulence kinetic energy 
cascades from larger to smaller scales. The cascade stops 
at a damping scale, Lp, where the turbulence energy 
is removed by dissipation. As usual, we assume that 
the turbulence is produced in a time much less than the 
Hubble time, T st i r <C )JH- m — here H m is the Hubble 
parameter at t- m — 0- 0] ■ an d therefore we ignore the ex- 
pansion of the Universe when studying the generation of 



GWs. In this case the GW equation of motion, in wave 
number space, is [l^j ] 



hij(k,t) + k 2 h l3 {k,t) = l&-KGU.ij(k,t). 



(1) 



Here G is the Newtonian gravitational constant, 
and hij(k) = J d 3 x e ik x hij(x) and hij(x) — 
J d 3 k e~ lk ' x hij(k)/ (2ir) 3 is the Fourier transform pair of 
the tensor metric perturbation which is defined hij = 5gij 
{ha = and hijk J — 0). We use natural units h = 1 = c, 
physical/proper wave numbers (not comoving ones), and 
an overdot denotes a derivative with respect to time t. 

Stochastic turbulent fluctuations generate stochastic 
GWs. Gaussian-distributed GWs may be characterized 
by the wave number-space two-point function 

{h^t)h lm {k! ,t)) 
= (27r) 3 5( 3 ) (k - k') [M ijlm H{k, t) + iAi jlm H(k, t)} . (2) 

Here H(k,t) and Tt(k,t) characterize the GW amplitude 
and polarization, 4Miji m ( k ) = PuPjm+PimPjl-PijPlm, 

and 8jAijl m (\t) = k.q(Pj m €Hq-\- Pi[Cj m q-\- Pim£jlq~\~ Pjl^imq) 

are tensors, and e^-; is the fully antisymmetric symbol. 
Choosing the coordinate system so that unit vector £3 
points in the GW propagation direction, using the usual 



circular polarization basis tensors e i 



(ei ± ie 2 ) 



(ei ± ie 2 )j/V2 ,13], and defining two states h + and h~ 
corresponding to right- and left-handed circularly polar- 
ized GWs, we have h%j = h + efj + h~e~j. The GW degree 
of circular polarization is given by |l4| 



V(k) 



(/i+*(k)/i+(k') -h-*(k)h-(k')) _ U{k) 
{h+*(k)h+(k>) + h-*{k)h-(k')) ~ H(k)' 



(3) 



Both H(k,t) and Tt(k,t) are obtained by solving 
Eq. (JTJ , and are related to (k, t) . For instance, an 
axisymmetric stochastic vector source (non-helical turbu- 
lent motion or any other non-helical vector fieldHnduces 
unpolarizcd GWs with \h+(k,t)\ = \h~(k,t)\ |1 H H; 
the presence of a helical source alters this situation. 

To compute the induced GW power spectrum one must 
have the source two-point function (U*Jk,t)Hi m (k' , t')). 
This is determined by the fluid velocity two-point func- 
tion. For stationary, isotropic and homogeneous flow the 



2 



velocity two-point function is 0, 0] 

«(k) Uj (k')) - (2TT) 3 5^(k-k')[P lJ P s (k)+te l /k l P H (k)}. 

(4) 

Here Ps(k) and -Pfj(fc) are the symmetric (related to the 
kinetic energy density per unit enthalpy of the fluid) and 
helical (related to the average kinetic helicitv (u-(Vx u))) 
parts of the velocity power spectrum 0, llfl . Causality 
requires P s (k) > \P H (k)\; see p. 161 of RefH0. 

However, in the case of interest here, the source of 
turbulence acts for only a short time r, possibly not 
exceeding the large-scale eddy turnover time (cor- 
responding to length scale Ls) — for self-consistency, 
however, we assume that the source is active over a time 
r = max(, stir , TS )S - resulting in a time-dependent 
velocity spectrum. To model the development of helical 
turbulence during the time interval (ii n) ifi = Un + t) we 
make several simplifying assumptions: 

(a) Turbulent fluid kinetic energy is present on all 
scales in the inertial range ks < k < ko- Here ks — 
2-k/Ls and ko = 2tt/Ld. We also assume that the en- 
ergy is injected into the turbulence continuously over a 
time r, rather than as an instantaneous impulse |5j, |6j . 

(b) Unequal time correlations are modeled as |5j 

K(k,tK(k',t')) = (2n) 3 6^Ck~k')[P lJ F s (k,t-t') 

+ ie ij ikiF H (.k,t-t% (5) 

where the t — t' dependence of the functions Fs and Fjj 
reflects the assumption of time translation invariance. 
Since energy is injected continuously, at t = t' € (t- m , £g), 
F s (k,0) = P s (k) and F H (k,0) = P H {k). 

(c) The decay of non-helical turbulence is deter- 
mined by a monotonically decreasing function D\ (t) and 
F s (k,t) = P s (k)Dt(t), p. 259 of Ref. Extending 
this assumption to the helical turbulence case we also 
model Fu(k,t) = Pj{(k)D 2 (t), where D 2 (t) is another 
monotonically decreasing function. Since in the consid- 
ered model most of the power is in the inertial range, for 
simplicity we discard power outside the inertial range by 
truncating P$ and Pjj at k < ks and k > kjj- 

(d) We model the power spectra by power laws, 
Ps(k) oc k ns and -Pff(fc) oc k UH . For non-helical hy- 
drodynamical turbulence the Kolmogorov spectrum has 
ns = —11/3. It has been speculated that in a mag- 
netized medium an Iroshnikov-Kraichnan spectrum with 
ns = —7/2 might develop instead. The presence of hy- 
drodynamical helicity makes the situation more complex. 
Two possibilities have been discussed. First, with a for- 
ward cascade of both energy and helicity (dominated by 
energy dissipation on small scales) one has spectral in- 
dices ns = —11/3 and nn = —14/3 (the helical Kol- 
mogorov (HK) spectrum), p. 243 of Ref. Second, if 
helicity transfer and small-scale helicity dissipation dom- 
inate [18|. ns = nn = —13/3 (the helicity transfer (HT) 
spectrum) 19]. The HK spectrum has been observed 



in the inertial ra nge of weakly helical turbulence (i.e., 
\Pti{k)\ *C Ps{k)) |20|. For strongly helical hydrodynam- 
ical turbulence the characteristic length scale of helicity 
dissipation is larger than the Kolmogorov energy dissipa- 
tion length scale 0, 0, 0] . Therefore the inertial range 
is taken to consist of two sub-intervals, both with power- 
law spectra. For smaller k the spectra are determined 
by helicity transfer and have ns = njj = —13/3, while 
for larger k turbulence becomes non-helical and the more 
common HK spectrum is realized. Since GW generation 
is mostly determined by the physics at small k 0- E| ■ it 
is fair to only use the HT spectrum in this case also. 

Based on these considerations we model Ps{k) = 
S Q k ns and P H {k) = A kg s ~ nH k n », where: (i) for the 
HK case 5 = 7r 2 C* fc £ 2 / 3 and A Q = Tt 2 C k S / '(e 1 / 3 ^) H3, 
implying A /So = S/(eks); and, (ii) for the HT case 
S = n 2 C s 6 2 / 3 and A Q = n 2 C a 5 2 / 3 T9.]. Here e and 6 
are the energy and mean helicity dissipation rates per 
unit enthalpy, and Cfc, C s , and C a are constants of order 
unity. 

Given a model of the turbulence, the turbulent source 
two-point function is 

(n*-(k,<)n im (k',i + 2/ )> 

= (2tt) 3 <5( 3 ) (k - k') [Miji m f{k, y) + iA l3 i m g{k, y)] , (6) 

where M.iji m and Aiji m are defined below Eq. J2|. The 
functions f(k, y) and g{k 1 y) that describe the symmetric 
and helical parts of the two-point function are 



f(k,y) 



^-/^/^^(k-Px-p.) 

x [(1 + 7 2 )(1 + p 2 )D 2 (y)P s ( Pl )Ps(P2)+ 

+ A 1 pD 2 (y)P H ( Pl )P H ( P2 )], (7) 

r , , (p + p) 2 Di(y)D2(y) f , 3 f , 3 
g(k,y) = j2g^g Jdpijdp 2 

x ^ 3 )(k- Pl -p 2 ) [(\+ 1 2 )l3P s {p 1 )P H {P2) 
+ {l + (i 2 ) 1 P H {pi)P S {P2)\ , (8) 

where 7 = k • pi and (3 = k • \) 2 . The helical source term 
g(k,y) vanishes for turbulence without helicity 

To determine H(k,t) and TC(k,t) we solve Eq. as- 
suming that there is no GW for times t < ii n , i.e., we 
choose as initial conditions /iy(k, t- m ) = = hij(\s.,t- m ). 
To compute the induced GW power spectrum we use the 
averaging technique described in detail in Sec. III.B of 
Ref. 0. The main points are: (i) the ^^(k - k') in 
Eqs. |(5J) and © ensure statistical isotropy of the GWs; 
(ii) the statistical average can be approximated by ei- 
ther a time or a space average (this is justified for lo- 
cally isotropic turbulence see Sec. 21.2 of Ref. |22|); (iii) 
we choose to time average since the Green function for 
Eq. Q and the source term n/ m (k, t) are time depen- 



3 



dent. These approximations result in Eq. (32) of Ref. 
(/i^(k,te)/i im (k / ,t fi )) ~ ^ 6 ^ k } T J dt cos(kt) 

x(m-(k,i 1 )n im (k',t 1 + t)>. (9) 



The two-point function on the r.h.s. of the integral is in- 
dependent of t\ and (h*j (k, ta)hi m (k' , ig)) is proportional 
to the source duration time r, as expected for locally 
isotropic turbulence, see p. 358 of Ref. |22l |. 

From Eqs. ® © we see that the symmetric 
H(k,t) and helical Ti.(k,t) parts of the GW two- 
point function in Eq. @ are integrals over y of 
cos(ky)D 2 (y), cos(fcy)Z^(y), and cos{ky)D 1 (y)D 2 (y). 
Taking into account that both Di(y) and D2(y) 
are positive monotonically-decreasing functions of y, 



ft* 1 dycos(ky)F a (p,y)F b {\k - p\,y) 



P a (p)Pb(\k 



p\)/(y/2k) (where a and b can be S or H). Integrating 
over angles, we find at t = tg, 

H(k) ~ Ajd PlPl J d P2 p 2 Q [(1 + 7 2 )(1 + #)P 5 (pi) 
x Ps(p2)+^p(3 p Ph(pi)Ph(p2)] 1 (10) 
- 2A / d Pl Pl J d P2 P2 Q [(1 +1 l){3 p P s (p 1 ) 

X Pff(p 2 ) + (l+^)7pi 3 H(Pl)i's(P2)] • (11) 

Here A = ar/(47r 2 /c 4 ) where a = ^2{p + P ) 2 (8ttG) 2 , 
7p = (k 2 +p 2 -p 2 )/(2k Pl ), P P = (k 2 +p 2 -p 2 )/(2kp 2 ), 
9 = d(pi +p 2 - k)6( Pl + k- P2 )8(p 2 +k- pi), and 9 
is the Heaviside step function which is zero (unity) for 
negative (positive) argument. 

For power-law Ps(k) oc k ns and Pn(k) oc k UH the in- 
tegrals in Eqs. (|10l) and illl|) can be done analytically 
but the results are complicated and do not edify. Instead 
we compute the degree of circular polarization, Eq. J2J, 
by evaluating the integrals numerically for different pa- 
rameter values. Results are shown in Fig. 1. 

Figure 1 and other numerical results show that for 
maximal helicity turbulence (when A = So) with equal 
spectral indices n^ = rig < —3, the polarization de- 
gree V{k) ~ 1 (upper solid line). For weaker helical 
turbulence (when Aq < So) with nn — ns < —3, 
V(k) — > CAq/Sq, where 1 < C(ns,nn) < 2 is a nu- 
merical factor that depends on the spectral indices. For 
HT turbulence with ns = nn = —13/3, C « 1.50, while 
for Iroshnikov-Kraichnan MHD turbulence (ns — nn = 
— 7/2), C w 1.39. Excluding the edges of the inertial 
range ks < k < kjj, an analysis of Eqs. i(TT)|) and ltTT|) 
shows that the main contribution to the integrals come 
from areas with p\ ~ ks, P2 ~ k and p\ ~ k, p 2 ~ ks- 
In this case (for arbitrary spectral indices ns, riH < —3) 
H(k), Ti(k) oc k ns ~ 3 kg S+3 , so this approximate analysis 
indicates that V(k) is not very sensitive to k. 




FIG. 1: GW polarization degree V(K,tf/), Eq. as a func- 
tion of scaled wave number K — k/ks relative to the large- 
scale wave number ks on which energy is pumped into the tur- 
bulence. This is evaluated at time ta, after the turbulence has 
switched off, and remains unchanged to the present epoch. It 
has been computed for a damping wave number ko — lOfes- 
Three pairs of curves are shown. Solid lines correspond to 
the amplitude ratio Ao /So — 1 (maximally helical turbulence 
,9]), dashed lines to Ao/So = 0.5, and dot-dashed lines are 
for Aq/So — 0.2. The upper line in each pair corresponds to 
HT turbulence with ns = nn = —13/3 |19ll21| and the lower 
line to HK turbulence with ns = —11/3 and nn ~ —14/3 
[l^|. Even for helical turbulence with Ao/So < 0.5, for large 
wave numbers k ~ kn, ns = nn = —13/3 is unlikely so the 
large K part of the lower dashed and dot-dashed HT curves 
are unrealistic. The large k ~ ko decay of the HK curves is 
a consequence of vanishing helicity transfer at large k 



In position space, we find 



9aT,S 2 fcg s+3 
"647r 4 (n s + 3) 



dk(l±V{k))k r 



(12) 

GW amplitudes are conventionally expressed as 
(ft«(x,tfi)^(x,tfl)) = 2j Q c °dlnf[\h+(f)\ 2 + \h~(f)\% 
see Eq. (11) of Ref. [23|, where the frequency / of a GW 
generated by an eddy of length L is / = 1 where tl is 
the eddy turnover time 0, UIH • This is inversely propor- 
tional to the cosmological scale factor, so the frequency 
/ today and / when the temperature was T m = 100 Tioo 
GeV are related by / = 1.65 x lir 5 Tioofl^o/ / H ia Hz 
[2^, where g- m = IOO5100 is the number of relativistic 
degrees of freedom at t- lri . Since we truncate turbulence 
power for L > Ls, the GW spectrum is non-zero only 
for / > fs, where 



/s = 1 , xl „- % /^(£)''V^\<"+ 5 >' 2 



\ns + 3| \vl \L S 



xT 100 H.-i Hz 



(13) 



is the frequency now that corresponds to the stirring 
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length Lg. Here we use the fact that for locally 
isotropic turbulence the energy dissipation rate e = 
2v f£ D dk k 4 Ps(k)/ir 2 (where v is the plasma viscosity, -~ 
p. 483 of Ref. 22]), is equal to the source power input, 
i.e., £ = 3Kpvac/(4pr) and k is the phase transition effi- 
ciency For HK turbulence with ng — — 11/3, Eq. I|13l) is 
Eq. (53) of Ref. 0. 

Using Eqs. I|12|) and (|13|l . and neglecting the weak k- 
dependence of the GW polarization degree V ', we find 
that h ± {f) oc /~ 11/4 for the HK case HQ, while for HT 
turbulence h ± (f) oc /~ 13 / 2 . We expect such a steeper 
dependence on frequency for helicity induced GWs, since 
the helicity transfer rate is more important on larger 
scales 0| . In both cases the amplitude of the GW spec- 
trum peaks at the stirring frequency fg. 

We close with a brief examination of the prospect 
of detecting such circularly polarized GWs. The GW 
energy density parameter for frequency /, Hgw(/) is 
given by (see Eq. (7) of Ref 0) n GW (f)h 2 = 5.9 x 
10 35 (|/i+(/)| 2 + \h-(f)\ 2 )(f/Rz) 2 , where h is the Hubble 
constant in units of 100 km sec -1 Mpc -1 . In our case, 



n GW (f)h 2 



1.05 x 10- 11 gJJ 3 



n s 



Ld 
Ls 



3(n s +5) 



3Kp vac Lc 



J2 



tH^ ' \ng + 3\ 

2(2n s +5)/(n s +5) 

(14) 



The stirring frequency fg and the GW spectrum are 
very sensitive to phase transition properties. If the phase 
transition is strongly first order, for the HK case fg ~ 
5 x 10 _3 Hz near the Laser Interferometer Space Antenna 
(LISA) frequency range, but the amplitude Ogw(/)^ 2 — 
10 -11 0, is below LISA sensitivity on this frequency 
0, ■ An additional limit of detectability is imposed by 
the dominating stochastic GWs signal from white dwarf 
binaries [24| . Thus it is unlikely that the GWs discussed 
here will be detected by currently planed GW detectors, 
but future detector configurations [2j| may well be able 
to. 
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